Bimetric varying speed of light theories and primordial fluctuations 

Joao Magueijo 

Theoretical Physics Group, Imperial College, London, SW7 2BZ 

(Dated: January 29, 2010) 

We exhibit a varying speed of light (VSL) theory that implements the recently proposed decaying 
speed of sound mechanism for generating density fluctuations. We avail ourselves of bimetric VSL 
theories, where the speed of gravity differs from that of light. We first show that a DBl-like type of 
A'-essence model has the necessary speed of sound profile to produce (near) scale-invariant fluctu- 
ations. We then examine the map between bimetric and if-essence models: typically the bi-scalar 
connecting the two metrics is a /("-essence field in one of them. Remarkably, the DB1 model is 
found to perturbatively represent the minimal bimetric model, where the bi-scalar is Klein-Gordon 
in the matter frame. But the full non-perturbative bimetric structure is even simpler: the bi-scalar 
dynamics should be simply driven by a cosmological constant in the matter frame, balanced by 
an opposite cosmological constant in the gravity frame. Thus the problem of structure formation 
receives an elegant and universal solution within bimetric VSL theories, which are known to also 
solve the flatness and entropy problems and evade a plethora of causality concerns. 

PACS numbers: 0000000 
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I. INTRODUCTION 

In a recent letter |l| we proposed a mechanism for pro- 
ducing scale- invariant density fluctuations of appropriate 
amplitude based on a decaying speed of sound c s . We em- 
phasized the mechanism's generality and how it could be 
implemented using a variety of methods. The examples 
of if-essence 0, Q and varying speed of light (VSL) [4|-|6| 
were given, and other possibilities, such as non-adiabatic 
hydro dynamical matter, were suggested. But we stressed 
the generality of the proposal and deliberately chose not 
to marry it to any specific model. It was nevertheless 
noted that important issues, such as the causality con- 
cerns raised by the presence of superluminal signals or 
the status of the other Big Bang problems (such as the 
flatness or entropy problems) did require the context of 
a specific model. In this paper we investigate specific 
models, with the focus heavily trained on varying speed 
of light theories. 

VSL theories may raise disturbing fundamental is- 
sues [T], |8|, [l9j , at worst breaking every symmetry one can 
think of. But this need not be the case; indeed they may 
be simply regarded as non-trivial realizations of the lo- 
cal Lorentz group, without introducing preferred frames 
and other "oddities" . Two implementations of this idea 
stand out: bimetric VSL theories and deformed special 
relativity. In the former the Lorentz group is realized 
by separate metrics for matter and for gravity, leading to 
different speeds for photons (and other massless "matter" 
particles) and gravitons [3, [lj| • In the latter the disper- 
sion relations are deformed from their usual quadratic 
form, so that the (group) speed of all massless particles 
becomes energy dependent [11| . To prevent the intro- 
duction of a preferred frame one then chooses a suitable 
non-linear representation of the Lorentz group [12j |. In 
two companion papers we examine the connection be- 
tween the varying speed of sound mechanism of [l| and 
these two types of VSL theories, concentrating on the 



bimetric approach here (see |14l4l6{ for connections with 
the latter). 

Bimetric theories are perhaps the most straightforward 
arena for realizing the varying c s mechanism for produc- 
ing scale-invariance. In |l| we exhibited an admittedly 
rather cumbersome X-essence theory with the required 
c s profile. Upon closer inspection, however, this model 
is found to be a limiting case of the Dirac-Born-Infeld 
(DBI) model (with the sign of one of its constants re- 
versed with respect to what is usual, so we should perhaps 
call it "anti-DBI"), as we prove in this paper. Further- 
more it has been shown that i\-essence models are re- 
lieved from causality paradoxes if interpreted as bimetric 
models [19|, [2Q| • In this paper we perform this exercise to 
discover that the DBI model universally associated with 
scale-invariance is precisely the minimal bimetric model, 
where the field linking the two metrics is provided with 
the simplest possible dynamics. 

The plan of this paper is as follows. We first review 
bimetric theories that fall under the VSL remit (Sec- 
tion|nJ), stressing a couple of details of technical relevance 
for this paper. We then recall the varying c s mechanism 
of [l| , highlighting the requirements for a viable theory of 
structure formation (Section llII[) . We examine the rather 
wooden i\-essence model previously proposed as "proof 
of concept" and show that this is in fact the limiting case 
of a DBI-like model (Section IIV|) . With this in mind we 
then consider the perturbative map between if-essence 
models and bimetric theories (Section [V|, showing that 
the DBI model considered maps onto a simple, minimal 
bimetric model. A few non-perturbative features of bi- 
metric theories are not captured by this map; but after 
identifying them we're able in Section IVII to write down 
the full non-perturbative bi-metric structure behind this 
model, which is simply a cosmological constant with re- 
spect to the matter metric. We conclude with a discus- 
sion of other possibilities motivated by this encouraging 
result. 



Throughout this paper we shall use a metric with sig- 
nature (H ); this is for consistency with the ma- 
jority of the literature quoted, but the reader is warned 
that some translating was required with regards to the 
DBI literature. We shall use the reduced Planck mass 
Mpi — l/(8irG) and set the present speed of light to 1. 



II. BIMETRIC VSL THEORIES 

Bimctric theories have been proposed as VSL theo- 
ries 0, [ly], solving the horizon, flatness, and dark matter 
problems. In the model proposed by Clayton and Moffat, 
for example, there are two metrics: an "Einstein" metric 
9fiv (which is used to construct the Einstein-Hilbert ac- 
tion and defines the "Einstein" frame, also called "VSL 
frame" in this context); and a "matter" metric g^y (used 
to construct the matter action via minimal coupling, and 
defining the "matter" or "constant c" frame). The left 
hand side of Einstein's equations is unmodified with re- 
gards to metric g pu ] but matter is minimally coupled to 

This set up is not unusual and has been used, for in- 
stance, in Brans-Dicke varying G theories [21(, where the 
two metrics define the Einstein and Jordan frames. What 
is peculiar to VSL bimetric theories is that the two met- 
rics are related by a disformal transformation, so that 
the light-cones associated with them don't coincide and 
the speed of "light" differs from that of "gravity" . TeVeS 
theories (an alternative to dark matter) can be seen as 
bimetric VSL theories [22]. 

The disformal transformation between the two metrics 
can be ruled by a choice of dynamics. In the simplest 
case one invokes a "bi-scalar" field 6 and sets: 



g^v = g^v + BO,, 



(i) 



Here B is chosen to have dimensions of M ~ 4 , so that <f> 
has dimensions of M. In the most general case B could be 
any function of 4> and X = -^d^d^cj)^ but in the minimal 
theory it is a constant. The action breaks into S = S g + 
Sm + Sa, with: 



On 




gR[g» 



gCmigfiv,® 



Matt 



(2) 



and S$ determining the bi-scalar dynamics. In the orig- 
inal formulation, due to Clayton and Moffat, <\> has a 
Klein-Gordon Lagrangian in the Einstein {g^ v ) frame, 
that is 



S,i. 



'''■'■v -A^d^d^-V) - C-ii 



But we could also make it Klein- Gordon in the matter 
{g^u) frame: 



Scf, — 




9 -g^d^d^ - V 



(4) 



or even consider hybrids, for which kinetic and potential 
terms live on different frames, e.g. 



S<:. 



d 4 x 



g l -r u d^3A- 



d 4 x, 



(5) 



The bi-scalar, after all, connects the two frames, so any of 
these dynamics is possible and can be considered minimal 
if B is constant. A subtlety, to be unveiled in Section IVTl 
reveals an even more "minimal" dynamics. 

We stress two little known peculiarities of these the- 
ories, which will be of technical relevance later. First 
the Einstein equations remain unmodified in the Einstein 
frame (G M „ = —T^/Mpj) only if we define the stress en- 
ergy tensor as 



T>"-= !_/ (S M + S#) 



(6) 



i.e. if we map the matter (and bi-scalar, if applicable) La- 
grangian into the Einstein frame before computing T^ v . 
The Bianchi identities are then promptly satisfied. But 
if we use the matter frame definition: 



rjip,L 



-{Sm + Sc/,) 
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(say, for theory (UJ)) then we have instead 



(7) 



QIH> = _ 



M" Pl 



' rpl_il/ 



(8) 



and consistency with the Bianchi identities is more intri- 
cate to prove Q • Cosmological perturbation calculations 
are most easily performed in the Einstein frame, so it is 
important to bear this in mind. 

A second peculiarity is that even though matter is min- 
imally coupled to g^v for both Lagrangian and equations 
of motion the same cannot be said of the bi-scalar. If the 
bi-scalar has a Klein-Gordon Lagrangian (either in the 
Einstein or matter frames) its field equations won't be 
Klein-Gordon, but receive a correction. This is because, 
upon performing variations in S, one finds 



Sg^u = Sg^ + 4:Bd {p (, 



(9) 



so that a term related to the stress-energy tensor appears 
in the field equation of cf>. If we take ([3]) this becomes 



V 2 + V'(0) 




D± M v A* v >■ 



= 



(10) 



This led [9J to the introduction of a "third" , derived met- 
ric, associated with the Klein-Gordon propagation of <j) 
(see their papers for details). If we take ((H) the field 
equation becomes: 



V 2 (f> + V'{<f>) 



BT^VpVrt 







(11) 



so that the "third metric" providing the field with a 
Klein-Gordon type of equation: 



r i "v /J .Vu<f>+v , {4>) = o 



T v = 9*" 



BT^ U 



(12) 



(13) 



This is not a proper space-time structure; for instance the 
covariant derivative V M used in (fT2"]) is still defined with 
respect to metric g^ v . But it renders the field equations 
in a Klein-Gordon format, something that will be of use 
later J44J . 

It is a notable fact that (in the absence of matter) the 
Lagrangian that does produce a Klein-Gordon equation 
in the g^ frame is a pure cosmological constant. This 
will be explained further in Section IVI1 but it should 
already be obvious from these expressions. 



III. THE DECAYING SPEED OF SOUND 
MECHANISM 

We now review the mechanism for generating scale- 
invariant fluctuations proposed in [l|. The central result 
is that for any (constant) equation of state w = p/p scale- 
invariance follows from a sound speed c s oc p, if the fluc- 
tuations originate from a vacuum state defined inside the 
(sound) horizon, according to the standard prescription. 

Let the speed of sound diverge with conformal time 
according to c s oc r)~ a (with a > 0; note that r\ is pos- 
itive and increases from zero). Whether we employ a 
hydrodynamical or a scalar field description, the density 
fluctuations are described by a modified harmonic oscil- 
lator equation. This can be written in terms of a variable 
v related to the curvature perturbation C by ( — —v/z, 
with z ex -S-. The equation for v is pi Il7l ITa] : 



c^k 

z 



v = 



(14) 



As with inflation this equation can be solved with Bessel 
functions, with solutions: 



(15) 



v = ^(3n(AJ„(l3c s k V ) + BJ-„(Pc s kr])) 
The order v is given by 

3(1 - w) 



f3 a 



2(l + 3tu) 



(16) 



and (3 = l/(a — 1) > 0. The boundary conditions, inside 
the (sound) horizon (kc s 7] S> 1) can be found from the 
WKB solution: 



ik I c s dr) 



-ij3c s kr) 



(17) 



These are satisfied if the integration constants A and B 
in (|15|l are k- independent numbers of order 1. The spec- 
trum left outside the horizon can now be found. Since 



c s rj is a decreasing function of time, the negative order 
solution is the growing mode, so that asymptotically we 
have: 



(18) 



frl 



{cskriY 



Scale-invariance of the curvature fluctuation (fc 3 £ 2 
const) therefore requires v — 3/2, i.e. 



a = «o = 6 



1 + w 
l + 3w 



(19) 



and if we rewrite c s in terms of the density p we conclude 
that this implies c s oc p for all w. The spectrum can also 
be made red or blue depending on whether a < c<q or 
a > cto, specifically 



ns - 1 = fi{a - a ) 



(20) 



Finally the fluctuations' amplitude can be worked out 
from formula (see Q): 



3^2 



k% 



(5 + 3w) 2 p 



1 



A I 



(21) 



Pl^-s 



giving further support for c s oc p as the root of scale- 
invariance. At low densities c s must be constant so scale- 
invariance must effectively require c s = Co(l + p/p*), 
where c s w cq at low-energy and p* is the density that 
triggers its divergence. Thus the normalization is set by 
the ratio of this scale to the Planck scale, viz: 



3a2 



fc J C 



(5 + 3u>) 2 



1 



w 



Ki 



10" 



(22) 



For the rest of this paper we shall use bimetric VSL the- 
ories to implement this mechanism, but for the sake of 
clarity it is important to dissociate our efforts from those 
previously made [23[ . 

In [23| a structure formation model was based on a bi- 
scalar <fi which is Klein-Gordon in the Einstein frame, i.e. 
satisfies © with V — 0. Therefore its speed of sound in 
the Einstein frame (where the fluctuations' calculation 
is done) is fixed, c s = 1. The fluctuations in follow 
Eqn. ((14)) with c s = 1 and decelerated expansion, and 
so their modes start off frozen- in (outside the horizon), 
oscillating only in the current epoch, when they enter the 
horizon. The model, thus, cannot rely on conventional 
methods for setting the initial conditions, i.e. by consid- 
ering (vacuum quantum) fluctuations inside the (sound) 
horizon and then following them as they leave the horizon 
and freeze-out. Quite the opposite happens: the <f> modes 
do not re-enter the horizon nowadays because they were 
never inside the horizon before. For the purpose of struc- 
ture formation the horizon problem for </> has not been 
solved. Thus the need to appeal to "non-conventional" 
initial conditions for super-horizon modes, as considered 
in El. 



Undoubtedly the horizon problem is solved for matter 
in these models. But matter is subdominant and irrel- 
evant for structure formation, as recognized in [23[. It 
may be that interactions between the bi-scalar and mat- 
ter causally connect super-horizon modes for the bi-scalar 
but these interactions are assumed to be sub-dominant 
(indeed matter is set to zero in [23j)- The main interac- 
tion of the bi-scalar is with gravity (according to (fT4")l 
with c s = 1), and fluctuation modes in the bi-scalar 
start frozen-in, not oscillating, stressing the presence of 
a "horizon problem" . 

With regards to "non-conventionality" , we do not feel 
qualified to throw the first stone. But we do believe it 
unnecessary to invoke new methods for setting initial con- 
ditions to obtain scale-invariance in bimetric models, as 
we demonstrate in this paper. 



subject to: 



p = -y=y/mMpi 



leads to solutions with p ex 1/i 2 and <j) = A/t. 
scaling solutions have a stable equation of state: 



w = \ - 



8M 



pi 



3 mA 



(28) 



These 



(29) 



but the model differs from the structure formation re- 
quirements in that c s — oo. In order to implement c s ex p 
a term of form L\ oc X n with n < 1/2 should be added 
to K. This doesn't affect the homogeneous solution and 
its equation of state at high energies, but c s is no longer 
infinite. Using Eqn. (f2l)f we find 



IV. THE (ANTI-)DBI MODEL AS A 
REALIZATION OF SCALE-INVARIANCE 

In [l| we already exhibited a proof of concept realiza- 
tion of the scenario just reviewed. Here we show that this 
realization is the limiting case of a Dirac-Born- Infold the- 
ory, with a crucial change in the sign of one constant with 
respect to what is usual in the literature. 

The model in [l[ is based on if -essence, a scalar field 
theory with non-trivial kinetic terms @, [2J-[26[ . The La- 
grangian has the form 



C = K{X) - V[<t>) 



(23) 



(with X = \d, 



2"MV 



!>), where K can be any function of X . 
Computing the stress energy tensor we find the pressure 
and density: 



P 
9 



K -V 
2XK X 



K + V 
whereas the speed of sound is found to be: 



K 



x 



Kjc + 2XK sXX 



(24) 

(25) 



(26) 



The cuscuton model [23, l28| is a good starting point for 
our construction. It is defined as a model with K ex 
\[X . It has a number of peculiarities: its speed of sound 
is infinite (since the denominator of (J26I) vanishes); its 
kinetic term K doesn't contribute to p but only to p, 
so that its energy density is fully due to the potential 
V if present. For a general potential the field can only 
take on discrete values if living on its own, but if other 
matter is present the cuscaton locks on to its dynamics. 
An interesting exception is a cuscuton with a suitably 
chosen mass term, which exhibits scaling solutions even 
without background matter [45j]. Specifically, 



C 



1 



2i2 



-m <f> 



(27) 



C Z S (X X5-" CX p 1 



(30) 



so that for n = —1/2 we have c s ex p, as required for 
scale-invariance. 

This argument may look contrived, but we repeated it 
here because it is "constructive" with regards to what is 
required from the point of view of producing fluctuations. 
But we now show that it can be directly motivated from 
the DBI action. As abundantly pointed out in the liter- 
ature [29h32| this model can be associated with "stringy 
physics" (an association which may or may not be rele- 
vant for this paper). DBI inflation has been studied in 
recent papers, where the model is tweaked so as to reduce 
the speed of sound during inflation. However the model 
could also be used to implement the varying speed of 
sound mechanism, dispensing with inflation altogether. 
The DBI model for a scalar field is based on the action: 



C 



1 



m 



y/l - 2f(<j>)X 



1 



m 



V{<j>) 



(31) 



where we note that we're using signature -\ (in 

agreement with all other literature used in this paper), 
as opposed to what has been used in the literature in 
this field. A positive / has the effect of limiting the 
curvature and reducing the speed of sound. We choose 
/ = — C where C is a positive constant, so as consider the 
minimal model with the opposite feature: the speed of 
sound increases at high energies, without any curvature 
capping. For X <C C _1 the theory has the appropriate 
limit C « X — V, so there's no question of introducing 
ghosts due to the signal of /. Usually X cannot grow 
above a given value (X = 1/2 /), but now it makes sense 
to explore the limit X ^> C -1 for which we obtain: 



£ = p 2 VX- 



/'" 



4y/X 



V(4>). 



(32) 



which is precisely the ad-hoc if-essence model we con- 
structed above (with p 2 = y/2/C), subject to a con- 
straint upon the coefficients of the terms in \/X and 
\j\[X. This finally lets us determine the value of p*, 



so important for the normalization of the spectrum, in 
terms of Lagrangian parameters: 



3 -A 

8 M 2 



F 



pi 



±01 

V2 w 



(33) 



In the latter identity we used ([2Uf to eliminate A in terms 
of w. The normalization implied by (|22p therefore re- 
quires that: 



(5 + 3w) 2 n A 



V2w(l 



w)M* l 



10" 



(34) 



relating field theory parameters and observables. 

Is there anything wrong with choosing the opposite 
sign for the DBI coupling constant? As explained above, 
there are no ghosts in the theory, provided a corre- 
sponding change in signs is introduced elsewhere in the 
Lagrangian (so that the low energy limit still comes 
out right). But the theory certainly doesn't come from 
string /brane theory, since the signature of the space fixes 
the sign of / to be positive (at least using a minimal 
adaptation of the argument). This is not a reason to 
discard it. The only problem with the theory seems to 
be exactly the feature we want to implement: a super- 
luminal speed of sound. This is a more general problem 
with some -ftf-essence models, which we now proceed to 
analyse. 



V. THE MAPPING BETWEEN A-ESSENCE 
AND BIMETRIC MODELS 

We now relate the two leitmotifs in this paper: bimet- 
ric VSL theories and the varying speed of sound mecha- 
nism as implemented by the DBI model, if-essence mod- 
els can be seen as theories with an emergent second met- 
ric [20] . Although the authors of [20] distance their work 
from full bimetric theories, in the limit of perturbativc 
fluctuations in a if-essence field the two perspectives are 
equivalent. So it's natural to enquire which emergent 
metric theory is associated with DBI models; and should 
we see such a metric as the signature of a bimetric theory 
(as we shall do in the next Section), which bimetric the- 
ory is associated with DBI models (or anti-DBI models; 
for the purpose of this discussion the sign of the coupling 
is irrelevant.) 

There is a good reason why if-essence and DBI models 
are usually only invoked when they reduce, rather than 
increase the speed of sound: the innate fear of faster 
than light propagation. Causality paradoxes, such as the 
anti-telephone, are expected to arise (see [T^, [20|; but 
also see [33| for an intrepid alternative). These can be 
circumvented by re-interpreting A"-essence as a bimetric 
theory. We review this mapping deriving the remarkable 
result that the DBI model capable of producing scale- 
invariance, as shown in Section IIVI is mapped precisely 
into the minimal bimetric theory discussed in Section [Til 



However we stress that the mapping discussed in this Sec- 
tion is perturbative, and an understanding of the full non- 
perturbative aspects of the bimetric theory will only be 
achieved in the next Section. 

The crux of the bimetric re-intcrprctation discussed 
in [20| is the remark that the field equation of a if-essence 
theory is equivalent to a Klein- Gordon theory subject to 
effective metric: 



M"" = K.xg^ + Kjix&Wt 



with inverse 



M K 



x 



,K. 



9^ 



xx 



K 



x 



d^4> .d v 



(35) 



(36) 



This metric is generally not conformal to the metric g^ v 
servicing gravity, immediately pointing us in the direc- 
tion of bimetric VSL theories. As long as 



K xx 



K 



x 



(37) 



(equivalent to requiring c 2 . > 0) the system is then hy- 
perbolic with respect to this metric and more generally 
the Cauchy problem is well defined, resolving causality 
paradoxes. So it's tempting to identify (1331) with a con- 
formal version of Jl} , with M M „ playing the role of matter 
metric <? M „ and action (H|). 

We stress, however, that this mapping with bimetric 
theories is imperfect. Firstly it induces a field equation 



M""V M V V <£ = 



(38) 



where the V operator is defined with respect to metric 
g M „, not M^ v . Secondly, a Klein-Gordon field equation 
for </) doesn't correspond to a Klein-Gordon Lagrangian in 
bimetric theories, as explained in Section [TTJ Finally the 
subtleties in defining the stress-energy tensor (see Sec- 
tion [n| are altogether ignored. Therefore the mapping 
just described can only be used gingerly. 

A stronger connection is established in [20j ]. consider- 
ing a background <f> — 4>q plus perturbations around it, 
a situation well suited to cosmology. Leray's theorem 
states that the perturbations satisfy 



M^DpDv S(j) = 
with a new metric (conformal to M^) given by 
w . K,xx, 



M^ = 



K 



x 



K 



x 



-vw< 



(39) 



(40) 



where all quantities are to be defined at <j) — cf>o. The 
covariant derivatives D^ v used in defining the Klein- 
Gordon equation are now defined with respect to metric 
M^ v and the theory may be derived from a Klein-Gordon 
Lagrangian fully built from M^ v . The mapping is there- 
fore rigorous and we should identify the inverse metric 



M^ = 



K 



x 



,K. 



9^ 



xx 



K 



x 



d^dvcj) 



(41) 



with the matter frame g^ v via a linearization of ([T]). The 
action for the perturbations 



1 



SS {2) = / d^xV^M-M^OMd, 



(42) 



is then understood from a linearization of Q or (JS|) via 
the proposed identifications (we shall ignore the potential 
for the time being, as only the kinetic term affects the 
mapping under study). 

With this mapping in mind we can now ask onto what 
-ftf-essence models one maps the minimal bimetric theory. 
In general any -KT-essence theory maps into a bimetric 
theory of form ([1]) but this could have a complicated form 
of B = B(<j),X). It is therefore interesting to ask which 
theories lead to a minimal bimetric theory. If we compare 
(HJ) with (|4"Tj) . ignoring the conformal factor, we get: 



_ 2 ^", xx K^xx 



K >x K x + 2XK. XX 



= B 



(43) 



This is equivalent to 

K t xx 
K.x 



B 



1 + 2BX 



(44) 



and if B is constant this integrates immediately into: 



1 



1 



K=-V1 + 2BX- u 



(45) 



which is the DBI model (notice that the potential does 
not enter the expression for the speed of sound and there- 
fore is left unconstrained). 

This result has a remarkable implication. If we con- 
sider a minimal bimetric theory where the bi-scalar is 
Klein-Gordon in the matter frame, then in a perturba- 
tive scheme we should identify g^ — Af M „. Then we find 
that in the Einstein frame, g^ Ul the theory maps into a 
DBI model, i.e. precisely the if-essence model that leads 
to scale-invariance. Comparing with Section IIVI we see 
that we have B = C, so that a positive constant B leads 
to a DBI model with / < 0, as required. This is a per- 
turbative result, but remarkable. It implies that scale- 
invariant fluctuations follow from the minimal bimetric 
model, with a constant positive B, with a minimal action 
defined in the matter frame. 

As a side remark notice that had we taken the origi- 
nal Clayton-Moffat theory, where the bi-scalar is Klein- 
Gordon in the Einstein (or VSL) frame we would have 
discovered that in the matter frame its dynamics maps 
into DBI model with / > 0. This can be proved by a 
straightforward adaptation of the above calculation, but 
the result can be intuitively understood. In this original 
model the speed of sound for the bi-scalar is the speed of 
gravity, which in the early universe is much smaller than 
the speed of light. So the speed of sound for the bi-scalar 
in the matter frame becomes negligible, just like for a 
DBI model with / > 0. 



VI. THE FULL NON-PERTURBATIVE 
STRUCTURE 

We are now prepared to identify the full non- 
perturbative bimetric structure behind the various mod- 
els presented in this paper. It will be simpler and more 
elegant, and perhaps we should have presented it up- 
front, bypassing the hardships of cosmology. We stress 
that the work in [20( is distinct from the bimetric picture 
to be presented in this section, and that there should be 
observational differences between the two, should there 
be further matter fields, or non-perturbative effects be 
under study. 

Let the dynamics of bi-scalar <fi be driven by a pure 
cosmological constant in the matter frame, that is: 



9 1 

'j,/, 



d 4 x 



ff(-2A) 



(46) 



Even if we don't bestow a kinetic term upon <f> (either 
in the matter or the Einstein frame) the field acquires 
dynamics, due to the subtlety explained at the end of 
Section |TTJ In the absence of K and V the first two terms 
in ([TTj) disappear, but since T^ v = Ag^ the field satisfies 



T^V^V, , 



A^Vp V* cj> = 



(47) 



that is, a Klein-Gordon equation in the matter frame. 
This is actually the only way to achieve a Klein-Gordon 
equation of motion in the matter frame, since a Klein- 
Gordon Lagrangian, Eqn. ((U), does not translate into a 
Klein-Gordon equation. 

In addition this theory has the required behaviour, for 
structure formation, to all orders in the Einstein frame. 
From (JT|) we have 



g = ,g(I + 2BX) 



(48) 



(see [3, |34j for a derivation) so that Si in the Einstein 
frame becomes 



2BX(-2A) 



(49) 



This is already all that is needed to obtain the correct 
high energy behaviour and therefore scale-invariant fluc- 
tuations. If additionally we want to ensure that the the- 
ory has the correct low energy limit we must choose 



A = - — 



1 
2B 



(50) 



and require that an exactly balancing cosmological con- 
stant is present in the Einstein frame. This results in the 
DBI action $T$ with / = -B (or C = B). As explained 
in Section IIVI we still need to add a mass term in the 
Einstein frame, with 



m 



2BM 2 pl 



(51) 



but this is in fact required by the Bianchi identities, if 
we want 4> to remain dynamical in all circumstances. 

To summarize all we need is a bimetric theory with 
two balancing cosmological constants, one in each frame: 



d 4 x 



1 
9 B 



d 4 x, 



(52) 



tuned to the same constant B which appears in the re- 
lation between the two metrics ([1]). Such a theory has 
a Klein- Gordon dynamics in the matter frame but be- 
comes DBI in the Einstein frame. The Bianchi identities 
then force it to have a mass term ensuring the appro- 
priate background scaling solutions at high energies. Its 
DBI behaviour in the Einstein frame induces a varying 
speed of sound with the correct profile to generate scale- 
invariant fluctuations with amplitude 



2 3 / 2 (5 + 3w) 2 



1 



w(l + w) BM% X 
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-10 



(53) 



(see Eq. IM|) . Thus, with the possible exception of the 
equation of state w, nothing remains undefined in this 
theory. 

It is possible that the reinterpretation of the anti-DBI 
action in the bimetric framework might shed some light 
on its known pathologies [41j if interpreted in a single 
frame (the theory only becomes anti-DBI in the gravity 
frame). This matter is beyond the scope of this paper. 



VII. DISCUSSION 

It is often stated that units can always be defined so 
that the speed of light is a constant. This is true in the 
same way that units can always be defined so that the 
Hubble "constant" is indeed a constant, the acceleration 
of gravity is the same everywhere as it is on Earth, the 
Universe is not expanding, etc, etc, etc. The argument 
is a perfect twin and likewise the reason it makes sense 
to talk about varying c is identical to why it makes sense 
to talk about the expansion of the Universe, Newton's 
laws, a varying G, etc, etc, etc: they are dimensional 
statements justified by the simplicity of the overall pic- 
ture. For example by choosing units where the universe is 
not expanding all the equations of physics would become 
ridiculously complicated. 

Nowhere is this more obvious, with regards to VSL, 
than in the present paper. We could always define units 
where c is constant: this is in fact the matter frame. But 
then the Einstein equations become a royal mess, with 



the perturbation equations looking even uglier; that the 
speed of gravitational waves is now changing is an omen 
of what would then happen to the formalism. The gravi- 
tational dynamics of the theory (background and pertur- 
bations) is crying out for the VSL or Einstein frame to 
be used, just as it would if we recast modern cosmology 
into units in which the Universe isn't expanding. This 
is far from a technicality. For example it could be ar- 
gued that the work in this paper is simply inflation [351 ] 
or ekpyrotic [361 ] in a different frame. This is obviously 
nonsense: if we erase the VSL feature from the theory 
by going into the matter frame Einstein's equations and 
all the perturbation equations are horrendously modified. 
This "dynamical" feature is not part of any inflation or 
ekpyrotic model, and no doubt the discussion of scale- 
invariance would be rather different if they were to be 
included Q. 

We have used -ftT-essence models, of the DBI persua- 
sion, as a stepping stone to the construction presented 
in this paper. But we believe that the bimetric VSL 
structure identified towards the end of the paper offers a 
superior description. To save ^-essence from the perils 
of the causality paradoxes the theory needs to be rein- 
terpreted as an emergent [201 ] or a bimetric theory. The 
construction obtained is then valid only perturbatively 
and coincides with ours to first order. But our bimetric 
construction is valid to all orders, and in any case has a 
simpler dynamics: two balancing cosmological constants, 
one in the matter frame the other in the Einstein frame, 
as discussed in Section IVTl 

In this paper we focused on initial conditions set by 
a vacuum state which is at first inside the horizon. We 
leave to further work establishin g th e connection with 
thermal initial conditions [l|, [37], [38] . This is because 
these require a phase transition in the speed of light, 
and are best realized in a theoretical framework diamet- 
rically opposed to the one described here (see, for ex- 
ample, [33, Eoj). The groundbreaking paper on VSL Q 
provides the best set up for implementing such a phase 
transition. In future work [42j we shall also present 
the predictions made by these models concerning non- 
Gaussianity and gravitational waves. After this paper 
was submitted new work appeared [431 ] where these issues 
are partly addressed. However the scope of [43| is differ- 
ent from our paper; notably inflation is still required. 

I'd like to thank N. Afshordi, R. Brandenberger, C. 
Contaldi, A. Tseytlin and A. Vikman for helpful com- 
ments. I'm very grateful to John Moffat for providing 
the inspiration for this paper and for a careful reading of 
the manuscript. 
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